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Abstract 

A self-consistent analysis of the currents in a tokamak yields the result that the parallel neutral- 
izing current is given simply by the sum of currents commonly called the bootstrap current and the 
Pfirsch-Schliiter current. An expression is given relating the parallel and perpendicular pressure 
gradient induced species flow velocities to the total current. 
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I. INTRODUCTION 



With a self-consistent analysis of the equilibrium currents in a tokamak, one finds that 
the parallel neutralizing current is given simply by the sum of what are commonly called the 
bootstrap current and the Pfirsch-Schliiter current. As Woods |l| correctly points out, the 
standard expression for the bootstrap current density 0, I^, Q], Jbs — {—ftrap/Bo){dp/dr) b, 
has a nonzero divergence, at odds with the usual property of source-less currents, V ■ J = 0. 
What Woods fails to mention is that the standard expression for the Pfirsch-Schliiter current 



density 



, Is], J PS — {—2ecos9/Bg){dp/dr) b, suffers the same fate. Together, however, they 



form the parallel neutralizing current density, Jii, which is without divergence. 



II. DERIVATION 



We begin with the usual equations of plasma equilibrium: 



Vp = JxB, /ioJ = VxB, V-B = V- J = B-Vp = J-Vp = 0, 



(1) 



where the total pressure is the sum of the species' pressures, p = J2s ^^Tg = neTe + {ni+nz)Ti, 
and J is the total current density. We know that a present-day tokamak has an ohmic current, 
Jn, and that we will be considering discharges with a neutral beam driven current, Jnbi- We 
work in a concentric circular flux-surface geometry {r,6,(j)), Figure 1, related to cylindrical 
coordinates {R,(j),Z) by i? = i?° + rcos^^ = + ecos9) , Z = —r sinO , = 0, where 
is taken along the direction of the plasma current. Various quantities also are defined 
in a coordinate system aligned to the magnetic field, {x,y,z) = (r, ±, ||), Figure 2, where 
z = B^/i?*^ = b = {0,bg,b(f,), and we assume toroidal symmetry, d/d(f) =^ 0. The magnetic 
field takes the form B = B°/(l + ecos^^) = {Q.Bq.B^) in this geometry. Some useful 
operators and identities [3| are: 
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(2) 



(3) 



Vx V/ = V-VxF = , V-(AxB) = B-Vx A-A-VxB , A-BxC = C-AxB. (4) 

We note that for nonzero J x B, = (dp/dr) f. The ohmic and NBI currents are solely 
in the toroidal direction, Jq + Jnbi = {Jn + Jnbi) 0, as is the measured total current 



3 = J (p, and satisfy V ■ (Jq + Jnbi) = 0. The sum of the driven and flow currents is 
J = Jq + Jnbi + J± + Jy, where J_l is the diamagnetic current and Jy is the neutrahzing 
current. The divergence of the diamagnetic current is given by 

V-J± = V- ^j^7^ ^ = (l/5^)Vp-VxB + B-Vpx Vx Vp = , (5) 

and we find that = p'B^/B"^ and = —p'Bg/B"^, for p' = dp/ dr. Turning our attention 
now to the parallel neutralizing current Jy, we find that in order to cancel the poloidal current 
coming from Jj^, we must have jj = —J^ = —p'B^/B^ . Then, by similar triangles Q], 
Figure 3, we have 



Bo ' 5° U 



4 = 'S'Je = ^(^) (l + .cos^).jf + Jf , (6) 



where we identify the first term as the bootstrap current and the second term as the Pfirsch- 
Schliiter current. That V ■ 3\\ equals zero is given by construction upon noticing that the 
divergence picks up only the poloidal component. So far, no mention has been made of 
"passing" versus "trapped" particles, but we note that only passing particles may participate 
in flowing currents. 



J = fi 



pass 



u" 



+ Jnbi , (7) 



for toroidal conductivity cr<^ = |&<^cr|| — bga±\ and pressure gradient induced species flows 
and ul. Thus, one may determine the profile given the total current, pressure, and 
magnetic field profiles. 



III. CONCLUSIONS 



That the bootstrap current and the Pfirsch-Schliiter current may combine to form a 
valid, source-less current density should not be a surprise, as the calculation by Ross[9] 
has indicated for some time. The magnitude of the diamagnetic and neutralizing currents' 
contribution to the total plasma current increases with the pressure gradient, suggesting 
that non-inductive current drive may suffice for sufficiently steep plasma pressure profiles. 
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FIG. 3: (Color online.) Magnetic Field and the Diamagnetic, Neutralizing, Driven, and Total 
Currents. 



6 



